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We study the full entanglement dynamics of two uniformly accelerated Unruh-DeWitt detectors 
with no direct interaction in between but each coupled to a common quantum field and moving 
back-to-back in the field vacuum. For two detectors initially prepared in a separable state our exact 
results show that quantum entanglement between the detectors can be created by the quantum field 
under some specific circumstances, though each detector never enters the other's light cone in this 
setup. In the weak coupling limit, this entanglement creation can occur only if the initial moment 
is placed early enough and the proper acceleration of the detectors is not too large or too small 
compared to the natural frequency of the detectors. Once entanglement is created it lasts only a 
finite duration, and always disappears at late times. Prior result by Reznik Q| derived using the 
time-dependent perturbation theory with extended integration domain is shown to be a limiting 
case of our exact solutions at some specific moment. In the strong coupling and high acceleration 
regime, vacuum fluctuations experienced by each detector locally always dominate over the cross 
correlations between the detectors, so entanglement between the detectors will never be generated. 



PACS numbers: 03.65.Ud, 03.65.Yz, 03.67.-a, 04.62.-fv 



I. INTRODUCTION 



Quantum entanglement, the uniquely quantum mechanical feature of a quantum system as Schrodinger emphasized 
[l|, is the distinguishing resource of quantum information processing (QIP). Its evolutionary behavior in time is thus 
of primary importance to the viability and functioning of quantum computing. Many QIP candidate systems involve 
qubits interacting with a quantum field, such as two level atoms (2LA) in a cavity. Because of the intrinsically 
OO 1 relativistic nature of a quantum field, in that physical information propagates at a finite speed, issues of causality are 
unavoidably imbued with QIP, even though the qubits in popular schemes are likely to move at nonrelativistic speed 
or remain stationary. It is of interest to ask whether quantum entanglement can transcend such limitations. This 
7-H ' underlies the frequently mentioned yet often-misconjured notion of "quantum nonlocality" such as usually associated 
<3\ with the famous EPR paper 0]. 

Following our two recent work discussing the behavior of quantum entanglement in a relativistic setting, one with 
, a detector (an object with internal degrees of freedom such as a harmonic oscillator or an atom) in relativistic motion 
! '"j ' [|[ and another focusing on the relativistic features of a quantum field mediating two inertial detectors [3| in this 
paper we study the conditions whereby quantum entanglement between two causally disconnected detectors (spacelike 
separated, outside of each other's light cone) can be created and if so how it evolves in time. By pushing to extreme 
conditions we can appreciate better this unique feature of quantum mechanics assessed in the more complete setting 
of relativistic quantum fields. 

The question in focus here is, whether quantum entanglement between two localized causally disconnected atoms 
can be created by the vacuum of a common mediating quantum field they both interact with, but not directly with 
each other. Refs. 041] affirm such a possibility whereas Refs. [1, ij| see no such evidence. These are not contradictory 
claims because the setups of the problem are not exactly the same. 

Using the time-dependent perturbation theory (TDPT) with extended integration domain, Reznik discovered that 
a pair of two-level atoms initially in their ground states will become entangled when they are uniformly accelerated 
back-to-back in the Minkowki vacuum of the field, even though in this setup the two atoms are causally disconnected 
Q. Later Massar and Spindel (MS) 0] considered an exactly solvable model with two Unruh-Raine-Sciama-Grove 
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(U-RSG) harmonic oscillators [10( moving in a similar fashion in (1+1) dimensions but both in equilibrium with 
the Minkowski vacuum. They discovered that entanglement can indeed be created, but only in a finite duration 
of Minkowski time after the moment that the distance between the detectors is the minimum. However, contrary 
to Reznik's perturbative results, MS found that such entanglement creation process does not occur in perturbative 
regime. In our assessment, since MS were looking at the late-time steady-state behavior (where "the detectors are 
in equilibrium with the Unruh heat bath" [3]), their result cannot be compared with Reznik's TDPT result, which 
corresponds to early-time behavior of the detectors or atoms. 

To resolve the difference and understand the discrepancy our present investigation adopts the physical system used 
in [i[ to a situation similar to that considered in 11 or 0.' Taking advantage of the existence of exact solutions to the 
model under study we can perform a thorough analysis and follow the system's evolution through their whole history 
which enables us to identify the conditions (in the motion of the two atoms), the approximations invoked (e.g., time 
dependent perturbation theory used in |5|) and the parameter ranges where quantum entanglement may be generated. 

This paper is organized as follows. In Sec. II we introduce the model and describe the setup in the problem. In 
Sec. Ill we derive explicitly the expression for the cross correlators, and discuss their evolutionary behavior. Using 
these correlators we examine the exact dynamics of quantum entanglement in different conditions in Sec. IV. We then 
give a comparison between the exact entanglement dynamics and the one from the reduced density matrix (RDM) 
of the truncated system. We conclude in Sec. V with some discussions. In Appendix A we give the expression for 
some elements of RDM from a first order time-dependent perturbation theory, and discuss some subtleties in the 
regularization and integration domain. In Appendix B we derive explicitly the RDM of the two truncated detectors 
in an eigen-energy representation. 

II. THE MODEL 

Consider two identical, localized but spacelike separated Unruh-DeWitt detectors, whose internal degrees of freedom 
Q A and Qb are coupled to a relativistic massless quantum scalar field as described in 1, under going uniform 
acceleration in opposite directions as described in |5[ . The action is given byEq.(l)in[3 , but now the trajectories of the 
detectors are chosen as z^ A = (a -1 sinh ar, a -1 cosh ar, 0, 0), and = (a -1 sinh ar, —a -1 cosh ar, 0, 0), parametrized 
by their proper times r and proper acceleration a. In this setup the detectors never enter into the light cone of each 
other, so they cannot exchange classical information and energy. 

Suppose the initial state of the combined system is a direct product of the Minkowski vacuum of the field and a 
separable state of the detectors in the form of a product of the Gaussian states with minimum uncertainty for each 
free detector, represented by the Wigner function 

p(Qa,Pa,Qb,Pb) = exp -\ (a 2 Q\ + a~ 2 P\ + fi 2 Q 2 B + /T 2 P|) , (1) 

where Pa and Pg are the conjugate momenta of Qa and Qb, respectively, and ln(a 2 /f2) and ln(/? 2 /0) are the squeeze 
parameters. Then the quantum state of the combined system is always Gaussian during the evolution since the action 
is quadratic. In this case the separability of the two detectors can be well defined by the covariance matrix of the 
detectors throughout their history. 

We study the entanglement dynamics by calculating the exact evolution of the two-point functions or correlators, as 
in To compare with the results in Q in the time-dependent perturbation theory regime, we will show the reduced 
density matrix (RDM) of the two detectors in an eigen-energy representation, then truncate the energy spectra of 
the detectors to include only their ground states and first excited states so that they look like two two-level systems. 
Using the explicit expressions for the elements of the RDM of the truncated system we can compare the criteria of 
separability in the dynamics of the exact and truncated systems. 

III. CROSS CORRELATORS 

When the detectors are in a Gaussian state the RDM obtained from integrating over the quantum field is fully 
determined by the two-point functions or correlators of the detectors. If we could obtain the time evolution of each 
correlator of the detectors, we have the full dynamics of this model. 

By virtue of the factorized initial state for the combined system and the linear interaction, each correlator splits 
into two parts as (• • •) = (• • ■) + (• • ■) . The a-part describes the evolution of the initial zero-point fluctuation in 
the detector, while the v-part accounts for the response to the vacuum fluctuations of the field. Since the quantum 
field effects such as retardation from one detector will never reach the other in the setup of this paper, no higher 
order correction from mutual influences is needed. So in this setup the expressions for the self correlators of a single 
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detector in [TTJ] and [|| are actually exact, and the v-part of the self correlators there can be directly applied here with 
( Qa(v) 2 > v = ( Qb(v) 2 ) y , ( Pa(v) 2 > v = ( Pb(v? > v , and ( Q A (r,), P A (ri) ) v = ( Qb(v),Pb(v) ) v - The a-part of the 
self correlators for Qa and Q b could be different if we take different values of a and f3 in the initial state ((T|) , but the 
calculation is still straightforward. The remaining task is to calculate the cross correlators between the two detectors. 

The a-part of the cross correlators always vanishes since the initial state ([1]) is factorizable for detectors Qa and 
Qb at the initial moment and no retarded mutual influence between them ever arises after the coupling is switched 
on. So the cross correlator ( Qa(ti), Qb{v') ) contains contributions solely from the v-part, or vacuum fluctuations of 
the field. Explicitly, it is obtained by performing the following two-dimensional integration, 



Qa(v),Qb(v') 



\2 

^Re 



ds 



(2) 



where Ao is the coupling constant between each detector and the field, 7 = Aq/87t, fl is the natural frequency of each 
detector, rj = t — tq, rf = t' — Tq are the durations of interaction from the initial moments To and Tq when the coupling 
with the field is switched on to the proper times r and r' of detectors Qa and Qb, respectively, and D + (z^(s), z' B (s')) 
is the positive frequency Wightman function of the massless scalar field, Eq. (|A11[) . In Appendix [S] we learned that 
one should take the value of the mathematical cutoff e in the expression (|AI1[) to be zero in calculating the cross 
correlators. Then the width of the \D + \ ridge on (s, s') plane (see Appendix IA 2[) is infinity in the A direction, and 
( Qa{t]),Qb{ti') ) can be written in the closed form, 
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where K- = (7 — ifl)/a, and Fr{x) = —2^1(1, 1 + K, 2 + K, x) X xj (I + K). If t , Tq — > —00 and 777, 777' -> 00 while 
ar, ar' are finite, only the first line of Q which is a function of r + t' will survive. This is the counterpart of the 
equilibrium result obtained by Massar and Spindel in the U-RSG model 0- 

Other cross correlators are obtained straightforwardly from ( Qa{i])tQb{i]') ) by proper-time derivatives, for ex- 
ample, ( Qa(v),Pb(v') ) = d T > ( Qa{v),Qb{v') ). and so on. 

Below we consider the case with r' = r and Tq = tq. In this case the time-slicing scheme is equivalent to Minkowski 
times. Given 
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as z — >• + , one can see that ( Qa(v),Qb(v') ) always vanishes as ar, a(r + To), and 777 — > 00. Thus the detectors 
must be separable at late times. But the transient behavior of the cross correlators could be more complicated. In 
particular, in the regime with 7 -C fl, a and 7 (r — To) -C 1 , but a(r — To) S> I, one can see that the cross correlators 
manifest the following multi-stage behaviors. 

In the cases with large — To > but still 7 |ro| <C I, ( Qa,Qb ) behaves differently in three stages respectively, as 
illustrated in FIG. HTleft): 

(i) When r < 0, the value of ( Qa, Qb ) is extremely small though exponentially growing (~ e 2aT ) if \a\ is not very 
small. 

(ii) When I < ar < — ar 0l the first line of fl3]) dominates, so ( Qa,Qb ) oscillates like 



Qa, Qi 



Tvye 



-2jt 



fl sinh 



7r irfl 

-It cos 2f2T H coth sin 2flr 

a a 



(5) 



where the r cos 2J7t term dominates at large t, when the amplitude of the oscillation grows almost linearly if 7T <C 1 . 
The amplitude of the oscillating ( Qa,Qb ) will reach the maximum at r ~ |to| if |to| < l/2 7 with the maximum 
amplitude about 
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FIG. 1: Evolution of ( Qa,Qb } for the detectors initially in their ground state, with 7 = 0.01, Q = 1.3, a = 2, % = 1, and 
starting at to = —60 (left) and to = (right). In the left plot one can see that when < r < —to = 60, ( Qa, Qb ) behaves 
according to ©, while after r > 60 it behaves according to (JS}. In the right plot one sees that for To = 0, ( Qa, Qb ) behaves 
like (J9j and its amplitude never grows. 



otherwise at r sa I/27 with the maximum amplitude about 
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(hi) After r > —To, the second and the third lines of ((3]) become important and cancel the tcos2£7t behavior, so 
one has 
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which oscillates with slowly decaying amplitude (« 2h^e~ 7T |to|/£! sinh(7rf2/a) if f2|ro| 3> 1). 

The underlying reason for the above three-stage behavior is similar to the one for explaining the behavior of pf[ 00 , 
which is discussed below (|A11|) . Later we will see that this three-stage profile of ( Qa,Qb ) will be present in the 
entanglement dynamics in the same regime. 

Two special cases are worthy of mention here. First, if 7T0 — > —00, ( Qa, Qb ) will never enter stage (hi). It will 
always behave like (J5J) at large positive r. Second, if tq — 0, ( Qa,Qb ) has no stages (i) and (ii). Rather, it starts 
with stage (in) and 
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which oscillates in proper time r with a very small and slowly decaying amplitude, as shown in FIG. [TJright). 

Note that the behavior of cross correlators is non-Markovian because it depends on the fiducial time, namely, the 
initial moment tq. So is the entanglement dynamics. 



IV. ENTANGLEMENT DYNAMICS 



A. full dynamics 

The degree of quantum entanglement between the two detectors in Gaussian states can be measured by the values 
of the logarithmic negativity E^f or the quantity S defined in Q. To demonstrate how the evolution of the cross 
correlators affect the entanglement dynamics, however, we calculate below c_ defined by 



Z± V^ 2 -4dctV 



1 1/2 
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with the covariance matrix V in (|B10|) and 

Z = dctvAA + detv_BB — 2detvAs- (11) 

c_ (c+) is the smaller (larger) eigenvalue in the symplectic spectrum of the partially transposed covariance matrix plus 
a symplectic matrix [J]. If c_ < h/2, the two detectors are entangled, otherwise separable. From c± one can easily 
determine the values for the quantity S = [c^_ — (7i /4)][c 2 _ — (h 2 /4)] we used before, and the logarithmic negativity 
by £V = max{0, — log 2 (2c_/?i)}, where the information in the range c_ > h/2 is beyond reach. 

Let us first consider the case with both detectors initially in their ground states, namely, a — /3 = y/Tl^, where 
2 + 7 2 is the renormalized natural frequency of the detectors. In FIG. [2] (left) we can see that the three- 
stage profile of ( Qa, Qb ) in FIG. [ljleft) emerges in the evolution of c_. A transient entanglement is created as the 
amplitude of the cross correlators grows, then decreases as the amplitude of the cross correlators decays, and totally 
disappears at a finite time. The created entanglement could remain in a duration much longer than the natural period 
of each detector. 

In the ultraweak coupling limit (7A1 <C a, f2) the feature of the cross correlators is even clearer in the entanglement 
dynamics. Indeed, in this limit we have (Til [T2| 

(Q A )~Q + O( 1 A ), (12) 



pD = {pr ^ 1 1 " 



^ 2 2+-^(Ai-ln^) +O( 7 A ), (13) 



with the constants Ao and Ai corresponding to the time scale of switching on the interaction and the time resolution 
of the detector, respectively, and 



Q= — 
20, 



e - 2 ^ + coth— (l-e" 2 ™) 
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(14) 



while ( Pa,Qa ) and ( Pb,Qb ) are 0(7) and negligible. When r is large, one can write the cross correlators as 
( Qa,Qb ) « xcos20r, {Q A ,Pb) = ( Pa,Qb > « -Oxsin20r, and ( P A ,P B ) w -^ 2 xcos2fir, where 

x = [t9{t) - (r + ro)e(r + To)] (15) 

a 

is the envelop of the oscillating cross correlators. Then one has 

c_ «fi(Q-|xl) + 0(7Ao,7Ai), (16) 

which is less than h/2 if the detectors are entangled. Now one can easily see how the profile of the cross correlators 
(~ \x\) enter in the entanglement dynamics. 

If l r o| 3> 1/7, the analysis is the simplest: When r > 0, Q has been in its late-time constant value and A term 
decays away. Then one has 

h 7rf2 A jh 



which, together with (J7J and (|12p. imply that there will be transient entanglement creation outside the light cone 
after r > if 



< 1 < 



ln(27fi/e 7 Ai) ~ Q ~ ln(e/2) 



10.238, (18) 



for the cases with Ai ^> | ln[7r/lti(27f2/e7Ai)]| in the ultraweak coupling limit. That is, to generate entanglement, 
the proper acceleration of the detectors or the Unruh temperature experienced by the detectors cannot be too small 
or too large, otherwise the self correlators will always dominate and the entanglement will never be created in this 
case. For the cases with the value of a/O satisfying (|18p. by substituting (|15l) into (|T7| . one can further estimate the 
moment of entanglement creation te ~ — Wq{— C)/^7 and the disentanglement time TdE ~ —W—i(— Q/2 r ), where 

, e-^/ a 7 /, , a\ . , nil 
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FIG. 2: Evolution of c_ of the detectors initially in their ground states (a = = y/TT^). (Left) The parameters in the left plot 
are the same as those in Fig. [TJleft) (i.e., to = —60, etc.) to facilitate direct comparison. One can see that the three-stage 
profile of ( Qa,Qb } emerges in the evolution of c_, and creates transient entanglement during 19 < r < 80, which is much 
longer than the natural period of the detectors. (Right) The right plot assumes the same conditions as the left plot except 
that the initial moment is to = —10. Now c_ is always greater than 0.5 because the cross correlators have not had sufficient 
time to become strong enough to create quantum entanglement. This shows clearly that the entanglement creation process is 
intrinsically non-Markovian in nature. 




FIG. 3: Evolution of C- of the detectors initially in their ground state (a = ft = \/Q, r ) in the ultraweak coupling limit, with 
7 = 10 -5 , Q, = 2.3, % = a = 1, Ao = Ai = 20, To = —2/7 (left) and to = —l/4q (right). Quantum entanglement is created 
(c- < H/2) after t > 0, but disappears (c- > h/2) at late times. 



Wk is the fc-th product log or Lambert W function, which is the inverse function of f(W) = VFexpVF. For example, 
for To — > —00 with other parameters the same as those in Fig. [3J one has T^wlx 10 3 and t^e ~ 2.8 x 10 5 . Note that 
in this setup the moment of entanglement creation te is always positive, when the two detectors are moving apart. 
While the 7A1 term in (jTTJ) is small compared to the Q — \x\ term, it can strongly affect the values of te and ts,e if 
g-Trfi/a is very small, i.e., f2(Q — Ixl) is very close to h/2. 

For the cases with a smaller |to| ~ 0(1/27) the situation is similar. Different values of To would give about the 
same te, if entanglement creation still happens, while the disentanglement time t^e and the minimum value of c_ 
(thus the upper limit of a/fl for entanglement creation) can be quite different but of the same order as those with 
I To I 3> 1/7- An example of entanglement creation in ultraweak coupling limit is given in Fig.|3l 

When a gets smaller than 7, the above approximation in ultraweak coupling limit fails [111 ]. The entanglement 
dynamics for the cases with both detectors at rest (a = 0) has been discussed in Ref. Q, while the separation d there 
should be taken as 2/ a ^ 1 in this setup. 

In the strong coupling limit, the self correlators always dominate over the cross correlators, due to the manifestation 
of the long-range autocorrelation in the field that each detector experienced locally in space. So quantum entanglement 
is never created if the coupling is sufficiently strong (see Fig. 0} . 

For the cases with the initial state as a direct product of two squeezed states of free detectors with minimum 
uncertainty, namely, (ce,/3) 7^ (V^r, the only difference is the a-part of the self correlators. Detectors in such 

cases will still be separable at late times, because all of the a-part of the correlators eventually decays away, and the 
late-time steady state is independent of the initial state of the detectors. In transient, the a-part of the self correlators 
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FIG. 4: Evolution of c_ of the detectors initially in their ground states in the strong coupling limit, with 7 = 0.1, = 1.3, 
% = a = 1, Ao = Ai = 20, and to — —60. Owing to the strong coupling the self correlators always dominate over the cross 
correlators , so quantum entanglement is never created. 




FIG. 5: Evolution of E with each detector initially in a squeezed state with minimum uncertainty. The parameters are the 
same as those in Fig. [T](left) except a and f3, which are indicated in the plots. In the dark spots of the above plots, E < and 
the detectors are entangled. One can see that the parameters in (a, 0) space far away from the ones for the ground state (at 
{a, 13) = (vHri \AV)> hi the third plot from the left in the upper row) tend to increase the value of E, namely, decrease the 
degree of entanglement. 



with the initial squeezed state oscillates in time about the one with the initial ground state: at some moment the 
former is greater than the latter, and at another moment, lesser. But the oscillations of ( ) a and ( Pj ) a are out 
of phase, so the overall effect of the a-part of the correlators is to increase the domination of the self correlators and 
decrease the degree of entanglement (associated with the value of £ or Ej^) if the values of (a,/3) are sufficiently 
far from (V^r, VTl^) (or the squeeze parameters are sufficiently large; See Fig. [5j. From our numerical results we 
observe, however, that it is still possible to enhance the tendency to entangle if one takes the values of (a,/3) very 
close to the ones for the initial ground state, while the enhancement is very tiny and the correction to the values of 
£ is usually of the next order in 7. 
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B. Dynamics with Truncated RDM 

At early times in the perturbative regime the leakage of amplitudes to higher excited states are negligibly small, so 
we expect the two detectors can be approximately seen as a two 2-level systems (2LS). For this kind of 2 x 2 system 
the two 2LS are entangled if the partial transpose of the RDM p„ Ang n , n , in (|B11[) has a negative eigenvalue, namely 
if one of the inequalities 

Pu),ioPoi,oi — Pn,ooPoo,ii — 0> (20) 

and 

Poo.ooPn.ii — Pio,oiPol,io — 0' (21) 

is violated. Reznik discovered [|| that the inequality (|2TJ)) for the RDM of the two 2LS will always be violated in 
TDPT with integration domain extended to R 2 (denoted as TDPTqo; Cf. (| A2|) - (j A5[) ) . according to which he claimed 
that quantum entanglement can be generated outside the light cone. 

In Appendix [B] we derive the truncated RDM for two UD detectors explicitly. Comparing this result with the full 
dynamics, we can make the following observations: 

1) Reznik's conclusion based on TDPTqo is not generic for all values of the initial moment To, the proper acceleration 
a, the natural frequency SI, or the coupling strength 7. Actually, the key result Eq.(18) in [5[ is correct only at the 
moment r = —tq in the TDPT regime that entanglement creation can occur. 

2) Beyond the TDPTqo regime, our result shows that once quantum entanglement is created, it only survives in a 
hnitc duration. The two detectors are always separable at r < and at late times. 

3) The inequality ([20|) implies that, if the quantum state of the detectors is separable, then from ([Blip one has 



< J AA ' J aa - J Aa J 



BB' ^AB rrA'B' _ ™* det(V + V 



16ft 10 



I> ( 22 ) 

Q 2 a)(Q 2 b)-(Qa,Qb' 2 



which is a product of the quantity E and a positive definite function of time. (Here Vb = (7j/2)diag(f2 -1 , fi, SI -1 , S7) 
is the covariance matrix for two free detectors in their ground states.) Thus the criterion of separability (|2T))) , though 
derived from the RDM of the truncated oscillator, is actually equivalent to S > 0, which is the sufficient and necessary 
criterion for Gaussian states in the full dynamics In other words, the RDM of the two oscillators each with a 
truncated spectrum, including only up to the first excited state, possesses the complete information pertaining to the 
separability of the two detectors in Gaussian states. 

We illustrate the time evolution of the exact Ipf^ool an d |Pio,io I lrom (|B11|) in Fig. |6l where the value of ool i s 
highly dependent on the behavior of the cross correlators (see Eqs. (|B22p - (lB26[) ): it grows as the amplitude of the cross 
correlators grows. At the same time the growing cross correlators tend to decrease the value of |pio,io I - Combining 
these two factors we find, in exactly the same duration that the detectors are entangled in the full dynamics, \pf x 00 | 
exceeds |pf 10 | and violates the inequality (|2"0|). 



V. DISCUSSION 



A. Summary of our findings 



We have studied the entanglement creation process of two uniformly accelerated UD detectors moving in opposite 
directions in the Minkowski vacuum of a massless scalar field. The two detectors are causally disconnected during the 
whole history and are far apart at late times. 

For two detectors initially in their ground states, entanglement creation does occur under some specific conditions: 
if the initial time is negative and not very close to zero, the coupling strength 7 between each detector and the field 
is not too large, and the ratio of the proper acceleration a to the natural frequency of the detectors f2 is at some 
moderate value. The moment of entanglement creation, if any, is always at positive r. Once quantum entanglement 
is created it can last for a lifetime much longer than the natural period of the detectors in some parameter range, 
while entanglement always disappears in a finite time. 

For two detectors each initially in a squeezed state, similar entanglement creation can also occur if the squeeze 
parameters in the initial state of both detectors are sufficiently small so the initial state is close enough to the direct 
product of the ground states. 



9 



0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 

-30 30 60 90 

FIG. 6: Evolution of |pn i( )o! (dotted curve) and | /°io, io 1 (solid curve) of the truncated RDM under the same condition as in 
Fig. [2](left). The duration the detectors are entangled determined by the condition |pn,ool > IPio.iol matches exactly with the 
result determined by c_ < h/2. 

Moreover, we find that the RDM of two oscillators each with a truncated spectrum up to the first excited state 
contains the complete information about the separability of the oscillators in Gaussian states. In Appendices A and 
B we see that in TDPTqo regime, where the integration domain has been extended to R 2 , the growing rates of the 
matrix elements of the truncated RDM can be estimated well, but one has to be careful in obtaining the ratio of 
different matrix elements. While TDPTqo results are Markovian, we find TDPT still keeps some of the non-Markovian 
features. 

The dependence of entanglement dynamics on the initial state, the fiducial time, and the non-Markovian features 
shown above are beyond the scope of Massar and Spindel's steady-state calculation In addition, there are some 
other interesting differences between our findings in (3+l)D UD detector theory and MS's in (1+1)D U-RSG model. 
First, in our analysis Ao and Ai corresponding to the ultraviolet cutoffs in (3+l)D UD detector theory explicitly 
enter the entanglement dynamics, while the infrared cutoff in (1+1)D U-RSG model does not affect the degree of 
entanglement in MS's result because they only considered the cases where their constant K, which is a product of the 
coupling strength and the very long duration of the interaction corresponding to the infrared divergence, is always 
much greater than all other parameters, such that the higher-order terms in the if-expansion of the logarithmic 
negativity can be neglected. Second, while in both cases quantum entanglement between the detectors can be created 
after the moment in Minkowski time that the distance between the detectors is a minimum, the lifetimes of the 
created entanglement in our case are usually much longer than the natural period of the detectors, in contrast to 
the short lifetimes of entanglement in MS's (1+1)D results. Finally, in our (3+l)D case entanglement creation is 
totally suppressed in the strong coupling regime and can manifest in the weak coupling limit, while in MS's result 
entanglement generation is a non-pertubative phenomena so there is no entanglement creation in their weak coupling 
regime. 

B. Quantum nonlocality, causality, and reference frames 

Entanglement creation between two causally disconnected objects as we have captured in this work may be viewed as 
a manifestation of "quantum nonlocality" in quantum physics. Quantum entanglement between two localized objects 
can be generated by allowing them to interact with a common quantum field even though they do not exchange any 
classical information. Foremost our results testify to the important fact that quantum nonlocality does not violate 
causality [l3j]. Notice that quantum entanglement can be recognized and put to use (as in QIP) only by those 
"spectators" who can access the information from both detectors. From the viewpoint of the separate detectors each 
can never find out the existence of the other from its own RDM, nor the quantum entanglement between them. There 
is no transmittal of physical information between them and causality (of information) is always respected. 

This fact also means that the existence of a "spectator" is essential when we refer to the dynamics of entanglement 
between the two localized parties outside of causal contact. The spectator is causally connected to (and thus can 
"see" ) both parties and can recognize the quantum entanglement between them through its own observation of both. 
Recall that in Ref. Q we found that entanglement dynamics for two quantum objects in relativistic motion depend 
on the choice of reference frames or coordinates (Minkowski or Rindler, for example). One may debate on which 
coordinate is "better" or more "objective" for the depiction of entanglement dynamics. Our results show that the 
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only physically meaningful one in a given physical setup for describing the entanglement dynamics of the system is 
that of the spectator. This is of course not unique, but there are well-established ways to relate what is measured by 
one spectator to another with considerations of time-slicing (see, e.g., Q) and transformation laws of reference frames 
in relativity theory. 



C. Quantum entanglement generated by collision 

The setup of our problem can also mimic the situation of two ions of the same charge in a head-on collision. 
Quantum entanglement is generated in the later stage of the collision when the two ions are moving apart, once the 
trajectories of the ions possess the right kind of symmetry to enhance the field correlations as those in the present 
problem (see Appendix IA 2} . There is no energy exchange between the two ions during the entanglement creation 
process of this kind. 

It has been shown that two atoms coupled with a common field vacuum in a cavity can get entangled in a "collision" 
[T^ . [l5| . These entangling processes could be interpreted as a consequence of virtual quanta exchange in a van der Waal 



potential, analogous to those with the coulomb interaction in electron-electron scattering. In 15[ the interaction time 
is much longer than the propagation time of photons across the spatial separation between the atoms, and energy 
exchange between two atoms is clearly present. These cases where entanglement generation occurring well inside 
the light cone is different from the present setup but closer to the situation we considered in where quantum 
entanglement is established mainly by retarded mutual influences. 

Note that particle (or atom, molecule) collisions in non-relativistic classical and quantum mechanics are usually 
described by an effective Hamiltonian with a direct, nonlocal, inter-particle interaction (e.g. a potential V(|ri — r2|), 
where rx,2 are the positions of the particles). Direct interaction always generates correlation between the particles, 
classical or quantum. The nonlocality of the interaction, on the other hand, is a consequence of coarse graining from 
a more fundamental local theory, resulting in an effective theory description. The particles in this non-relativistic 
regime cannot resolve the time scale for information propagating back and forth between them. Thus the generation 
of correlation by this kind of nonlocal direct interaction happens in a time scale during which two particles have 
causal contact for a long time. One cannot tell whether correlations can be created between two causally disconnected 
particles in this kind of effective theories. Moreover, with this kind of nonlocal direct interaction, if one does not 
further introduce a spatial range of interaction associated with a coarse-graining in time, even classical correlation 
will be nonlocal and cannot be described by any local hidden- variable theory. 
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Appendix A: Matrix elements in time-dependent perturbation theory 

In TDPT, one expands the wave function of the two detectors in terms of the energy eigenstates <p m (QA) and 
Vn(Qs) of the free detectors as 

1>{t)= ]T C m/n/ , mn e-^' + ^ M ^ + ^ T ^ m ,(Q A )^ n ,(Q B ) 7 (Al) 

then calculates the first few factors in the perturbative series of C m > n >, mn assuming a small coupling constant Ao, 

namely, C m > n ^ mn = J2n ^m'li' mn- Here oj n = E n /h = Q[n + (1/2)]. For the initial state JT]), to lowest order in 
Aq , the elements of the RDM of the detectors are 



pfo.io * 2MI [ ds [ ds '^ iU{s ~ S ' )D+ ^A{s),z v A {s')), (A2) 



ds ds'e- l0(s - s ' ) J D + (z^( S ),z^(s')), (A3) 



p !(),oi ~ / as J ds ' L> ' u .v- s) - : i>- 



TO •> TO 



Pn.oo « -^7^ / ds / ds'^^D+(z^(s),zUs')), (A4) 
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Pii.ii 



rfsds'dsds'e^^'- 5 - 5 ') [D+ (z A {s),z B (s'))D+ («£(*),«£(*')) + 



2/ift, 

i?+ (z B (s'),z B (s'))+D+ (z A (s),z B (s>))D+ (z B (s'),z»(s))} , (A5) 

etc., where D + is the positive frequency Wightman function of the massless scalar field in Minkowski vacuum, given 
by 



D+(z?,z»,) = (O m |0(^)0(4)|Om) 



M 3 fc _ we iaJ ( z 9_ z o ) +4k .( Zj _ z ) 



(2tt) 3 2w 

^/4tt 2 



,2 (A6) 

2 _ ,0 _ ,0 



with j, j' — A, B and the frequency uj = |k| for the massless scalar field, e is a small positive real number serving as 
a regulator at high frequency in k-integration. 

1. The regulators 

In obtaining pf 010 in (IA2[) . one could insert the trajectory of the detector A into (|A6|) . which gives the Wightman 
function as 

D+(z%(s), z v A (s')) = _^ sinhfA(sinh l / ^_ i£acoshaT) +£2 > (A7) 

with T = (s + s')/2 and A = s — s'. Conventional wisdom says that the value of e is extremely small, so the expression 
of the Wightman function could be replaced by [l6[ 

D'+(z^s),zl(s')) = 4 ^ 4 ^ 2 



^sinh^f(A-ze') 
^/4tt 2 



- 4f sinh § A (sinh § A - ie'a cosh § A) + e' 2 cosh a A ' 



(A8) 



which is independent of T. Indeed, when aT is small, the original D + and the modified D /+ have similar positions of 
poles and similar behavior around the poles in the complex A plane. 

However, as shown in Fig. only the modified D l+ will give a linear-growing phase in transient like 

Puvo = 27^/(e 2 ^ /a -l), (A9) 

where the growing rate PicfioA? can be obtained by extending the integration domain of (s,s') to R 2 , namely, from 
TDPTqo. In contrast, the original D + leads to a totally different behavior of /?{g 10 : It is linearly decreasing in the 
beginning and saturates after r > — a -1 Inc. 

Learning from what we did in taking the coincidence limit of the self correlators [11]] , we find that, if we put the 
small shifts eo and ej at the upper and lower bounds of the s and s' integration domain, such as 



A 



2 fT+ei 



P^w « ^ / ds ds'e-^°-^D+ , (A10) 

™' -'to Jro+eo 



and take the e in D + to zero limit, then the exotic behavior of D + will be altered and we recover the results obtained 
by using the modified D' + with e' corresponding to eo when eo = t\. So we interpret e as the mathematical cutoff, 
which should be taken the zero limit at some point of calculation, while the eo, t\ in (|A10[) and the e' in (|A8I) are 
physical cutoffs corresponding to the time resolution of the detectors and related to the constants Ao and Ai in the 
expressions for the self correlators [ll[ . The replacement from (|A7[) to (|A8[) actually changes the interpretation of the 
regulator. 
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FIG. 7: Early-time evolutions of the pf Q 10 in TDPT with the original positive frequency Wight man function, D + (z^(s), z^(s')) 
in Eq. (|A7[) (dashed curve) with e = e -14 /Q. is the Euler's constant), the pio.io with the modified one, D' + (z^(s), Za(s')) 
in Eq. (fA"8)) (solid curve) with e' = e~ 14_7c /f2, and the pg, jl0 in (TATO)) with the original D + (z%(s), z^(s')) and e = e~ 30 /ft, 
together with the physical cutoffs eo = ei = e _14_7c /Sl 3> e (dots). The left and the right plots are the same results but shown 
in different scales to stress that the early-time behavior of the dotted curve really matches with the climbing solid curve. Here 
7 = 10~ 5 and = 2.3. The slope of the dotted curve in its linear growing phase agrees quite well with 2y/ ( e 2nQ / a _ 1) m (|A9[) . 



2. Calculating |/9u,ool 

Usually in obtaining p^ 10 , TDPT is good from r\ = r — r » A±/Q up to r\ ~ 0(1/7) [ll| with 7 = Aq/87t (here 
the masses of both harmonic oscillators in the detectors have been set to mo = 1 and the proper accelerations of the 
detectors are not extremely large) . This is a large duration if 7 is small, which is consistent with the assumption of 
TDPT, so it is commonly argued that extending the domain of s integration from (to, t) to (—00, 00) in this duration 
would not change the results too much, that is, the TDPTqo results should be very close to the TDPT results. 

In Ref-H, Reznik compared the absolute values of |pf 10 ol an d |pfo,iol m TDPToo and discovered that the former 
is always greater than the latter, then drew the conclusion that the two detectors are always entangled. Nevertheless, 
with the extended integration domain, both \pf\ 00 | and |pfo 10 | are infinities, which are meaningless as elements of 
any density matrix. Usually one divides \pi$ 10 1 by the infinite duration of interaction and interprets the result as 
the transition probability per unit time in the linear-growing phase in TDPT regime [16j |. But the choice of proper 
duration for 00 | and the interpretation for the result could be tricky (see Appendix [Bl for more details). Further, 
our numerical results show that if we calculate \pf 1 00 | with a finite integration domain and a nonzero regulator e, the 
value of (|A4I) after integration can be quite different from those obtained in Ref.[5|. 

To calculate 00 , one substitutes the trajectories z^ A and z^ into the positive frequency Wightman function and 
obtain 

D+(z^s),zUs')) = — —, - ^ 4 f. r-r-TJXYT r ( AU ) 

^ cosh aT (cosh aT + tea smh | A) + e 2 

Now D + (z^(s), Zg(s')) is regular everywhere in (s, s') plane. For small finite e, since cosh aT > 1, the e 2 term in the 
denominator can always be neglected, but the e term cannot. If e sinh aA/2 cosh aT, then D + w a 2 /167r 2 cosh 2 aT, 
which is a function of T only and peaks around T = 0. For A sufficiently large such that ea sinh aA/2 > 1, however, 
\D + \ will be suppressed. So \D + (z^(s), Zg(s'))\ looks like a ridge at T w and spread with the width 4ri in A 
direction, where T\ = (1/a) sinh _1 (l/ae) w (1/a) ln(2/ae) (see Fig. [S]). 

Suppose the coupling is switched on at the initial moment r < 0. From the landscape of D + (z^(s), Zg(s')) shown 
in Fig. [51 the evolution of \pf x 00 | will have three stages: 

(i) When tq < t < 0, D + (z^(s), Zg(s')) is small all over the s-integration domain (within the square with two sides 
in long-dashed lines in Fig. [8] (right)), so the change in \pf[ 00 | is tiny. 

(ii) At about t = 0, the edge of the s-integration domain [within the square with all sides in solid lines in Fig. [8] 
(right)] touches the ridge of \D + (z^(s), Zg(s'))\ around T = 0, which gives obvious contribution to the result. As t 
grows, more and more portion of the ridge are included in the integration domain. Since the height of the ridge is 
almost constant in the A direction, |pf\ 00 1 grows linearly in r. 

(iii) The linear growth terminates at about t — t\. Then the integration domain has covered almost all the ridge 
so that \D + (z^(s), Zg(s'))\ does not add obvious contribution after t = t\ and \p^ x ol saturates at some constant. 
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FIG. 8: (Left) An example of \D + (z^(s), Zg(s'))\ in Eq. (|A1 1[) . Here e = e _s and a = 1. One can see that the correlation 
of the massless scalar field at z^(s) and Zg(s') is enhanced around s' = — s. Such an enhancement will be suppressed if the 
trajectories of the detectors are off the z^(s) and z B (s') in this paper or the scalar field is massive. (Right) Contour plot of the 
same \D + (z^(s), z B (s'))\. The long-dashed lines, solid lines, short-dashed lines, and dotted lines are borders of the integration 
domain at r = —10, 0, l/o sinh -1 e _1 (~ 8.69), and 15, respectively. 



From the argument in Sec. IA 11 one should take e — > so T\ and the width of the ridge in the A direction is virtually 
infinite. Then stage (ii) will be terminated around r w |to| because only a portion of the ridge with |A| % 2|rrj| is 
covered in the integration domain. This could be enough to make the value of \pf{ Q0 \ exceeds \pf- Q 10 | and the two 
detectors get entangled, if |to| is sufficiently large. On the other hand, for To > 0, the linear growth of stage (ii) never 
occurs, and |pf 10 ol i s always less than \p^ 10 |. 

Thus one can see that the value of \pf{ 00 | as well as the ratio \pf{ ol/|pfo 10I have complicated evolution in time 
even in TDPT. Their behaviors are quite different from the constant ones obtained in Ref.[5(. 



Appendix B: Truncated reduced density matrix of the two detectors 



Generalize the method applied in [11] for one single detector, the elements of the RDM in eigen-energy representation 
can be read off by comparing the coefficients of the s™ 3 terms in both sides of the following equation, 



E 



i \ nA 



GdetG 



lexp]>> [KiiG-^Kj 



5 lJ 



where i,j = A, A', B, B', K A = K A > = K B = K B , = K = y/TQh, Q is defined by 



5 = 4 



Qa)(Q 



and G is parametrized in 
/ 



G = 



K 2 



a A + 2 
a A - 
a x + + ib x + 
\ a X -+ib X - 



with the factors 



a A ± = 



■2Q 



a A - a x+ + ib x+ 

K 2 , 

a A+ - ib A + -f- a X - - ib X - 

K 2 

a X - - ib X - a B + +ib B + 
a x+ - ib x+ a B _ 



V u ± A (V- 1 ) 22 detV 



a B - 



a X - + ib X - 
a x+ - ib X - 
a B - 
ibf 



Ki 



(Bl) 



(B2) 



(B3) 



(B4) 
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a B ± = 



1 

2G 

1 

~2G 



V ™± 4 . 1)44 
ft 2 V 7 

yi3 ± ^(y-i) 24 det y 



fcx± = t^[((P.4,Qa)(Qa,Qb)-(P.4,Qb 



— ((P a ,Qb)(Qa,Qb)~(Pa 1 Qa)(Q 2 b)), 
^{(Pb,Qa)(Qa,Qb)-(P B ,Qb)(Q 2 a )), 



ng 

(( )(Qa,Qb)-( Pb,Qa ) ( Ql ))] • 



»± 



(B5) 
(B6) 
(B7) 
(B8) 

(B9) 



These a and b factors are obtained by solving the coefficients G y in the Gaussian RDM of the detectors p R 
exp — j QiG lJ Qj in terms of the two-point correlators. Here V lJ are the elements of the covariance matrix 



V = 



(BIO) 



w 1 vj 



in which the elements of the 2x2 matrices vy are symmetrized two-point correlators v„ 

( {n^n v n +n v n n^) ) /a with = (Q m {t),p m (t)), n,v = 1,2 and m,n = a,b. 

Up to the first excited states of both detectors, the truncated RDM of the two detectors in eigen-energy represen- 
tation reads 
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2 rrAB 





gK 2 J BB ' gK 2 J A ' B 
gK 2 J AB ' gK 2 J AA ' 














j-AB J-A'B' jAA' jBB' 



\ 



jAB' jA' B 



(Bll) 



with riA, fiB, n*, = 0, 1, = G , and 



V^detG 



(B12) 



Here we use the basis (riA, "s) = {00, 01, 10, 11} for p R . 

What TDPTqo could describe well is the transient behavior during ft -1 <C r < |r | <C I/27 while ft -1 <C r — t <C 
I/7 with to < 0, when the amplitudes of the oscillating cross correlators are linearly increasing, while the magnitudes 
of the self correlators are still in transient and growing linearly, too. In this regime, ( Qa, Qb) can be approximated 
by ([5]), which also yields 



(Pa,Qb) 
(Pa,Pb > 



Qa,Pb 



sinh — 



ftsinh^ 

a 

2ftr cos 2ftr — 



, 7rft 7rft . 

2ftr sin 2ftr + ( coth 1 ) cos 2ftr 

a a 



7rft 7rft 

— coth 

a a 



1 sin 2ftr 



up to 0(7). Together with the 0(7) results of the self correlators in TDPT regime from [TT[ with mo = 1, 



?A 



Qa,Pa) = (Qb,Pb) « ^ 



p2 



p2 



ft 

m 

Tr, 

o" 

2ft 
Tift 
"2" 







^coth 








^coth 











2 
Trft 1 



sin 2ftr] 



2ft 



n 7 ft 



Trft 



coth 1 ) 77 H — — [Ai - In a + (A - In a) cos 2 ftr?l - 

' 7rft L 



sin 2ft?7 
2ft 



(B13) 
(B14) 



(B15) 
(B16) 

(B17) 
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one finds that, for the detectors initially in their ground states, 



a B + £ 


Q 

* 2h 


fi, [ 7T 


a B - f= 




2Vtri 


ft 




- 1 


ax+ s 




7 


2Vt 


Tisinh^l 

a 


b x + s 




7 


-2 


ftsinh^ 





[ — In a 

1 /A 

+ -(A - 

7T 

2f2r cos 2f2r 



(Ao — In a) cos 2Qry] 
+ Ai - 2 In a) 



7rf2 7r£l 

1 coth — 

a a 



1 



-2fir sin20r- 



7rf2 7rf2 

1 coth 

a a 



• sin 2f2?7 



sin2i7r 



cos 2f2r 



up to 0(7), so that 



Ph. 00 



P01.01 



-9^ 



(^ + a A +) 2 



(B18) 
(B19) 
(B20) 
(B21) 

(B22) 
(B23) 



Now the growing rates of \pf[ 



00 



and \pq[ oil at large r agree well with those estimated by TDPToo, where both r\ 
and t are replaced by infinities. 

Nevertheless, one has to be careful in calculating the ratio of different matrix elements in TDPToo. Indeed, at 
nr > A 0) Ai, |pft. 

oil i s approximately proportional to rj — r — r while |pf 10 ol i s approximately proportional to r, 



so that \p R l ool/IPio 10 1 ~ e x 2r/?7, where r and 77 do not simply cancel each other. |pjj 00 1 / 1 Pio 10I i s actually 
time-varying, and Eq.(18) in Q is correct only at the moment 77 = 2t, or t = —To, which is almost at the end of the 
TDPToo regime. 

Note that, in weak coupling limit, 



a x + 
bx+ 
a B - 



Qa, Qi 



8(Q 2 A )(Q 2 B 



1 

1 



Pa,Qi 

(Ql) 



Qa-Pb ) 

(Ql) 



8(Ql 



2hr 



p 2 



(B24) 
(B25) 
(B26) 



so the behavior of o in (|B22[) is highly dependent on the behavior of the cross correlators, while p^i 01 in (IB23[) 
is not. 

Note also that (|B11|) in ultraweak coupling limit is numerically consistent with the TDPT results in Appendix [A] 
in the time interval fi -1 -C r — tq -C I/7 for < —To -C I/27. Here r can be negative or greater than |to|. From the 
three-stage behavior of the |pf x oo I in Appendix [3] which depends on the fiducial time t , one can see that TDPT 



still keeps some of the non-Markovian features, though the damping behavior (<■ 



,-2 7 T 



) has been lost in the coupling 



constant expansion (cf. the integrands of (|A2|) - (|A4|) with ([2])). On the other hand, the validity range of TDPToo is 
more restricted. TDPToo is valid only when both \p xl 00 | and \pf{ 00 \ are linearly growing, namely, only in the section 
fl^ 1 < r < |r | < I/27 of the time interval fl" 1 < r - r < I/7 for TDPT. TDPToo results are Markovian because 
even the memory of the initial moment is lost after the integration domain is extended. 
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